Gap function symmetry and spin dynamics in electron-doped cuprate superconductor 
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An antiferromagnetic (AF) spin fluctuation induced pairing model is proposed for the electron- 
doped cuprate superconductors. It suggests that, similar to the hole-doped side, the superconducting 
gap function is monotonic d^2_y2-wa.ve and explains why the observed gap function has a nonmono- 
tonic d^2_y2-wave behavior when an AF order is taken into account. Dynamical spin susceptibility 
is calculated and shown to be in good agreement with the experiment. This gives a strong support 
to the proposed model. 
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I. INTRODUCTION 



Pairing symmetry is an important issue towards un- 
derstanding the mechanism of superconductivity. For 
hole-doped cuprate high- Tc superconductor, it is gen- 
erally accepted that the pairing symmetry is (1^2 _y2- 
wave [ij. On the other hand, although no consensus 
has been reached yet, more and more recent experiments 
have pointed out that the order parameter of electron- 
doped cuprate superconductors is also likely to have a 
(ia;2_j,2-wave pairing symmetry [3, 0, 0, S @]- Inter- 
estingly however ^ngle resolved photoemission measure- 
ment (ARPES) [1| and Raman scattering 0, Q suggest 
a nonmonotonic dx2-y2-'wave order parameter with max- 
ima close to the nodes (diagonals) rather than to the Bril- 
louin zone (BZ) boundary. Understanding the origin of 
the nonmonotonic da;2_j^2-wave order parameter becomes 
an important issue. Yoshimura and Hirashima performed 
a strong-coupling one-band calculation and claimed that 
the nomonotonic feature comes from a strong AF spin 
fluctuation In contrast to hole-doped cases, they [?} 
found that the hot spots, the intersections of the mag- 
netic BZ boundary and the Fermi surface (FS), are lo- 
cated near the diagonals of the BZ. Alternatively it was 
also argued that the nonmonotonic feature of the order 
parameter is the outcome of the coexistence of the super- 
conducting (SC) and the AF orders When AF 
order coexists with the SC order, the resulting quasipar- 
ticle (QP) excitation can be gapped by both orders and 
behave to be nonmonotonic dx2-y2-'wave, even though 
the SC gap itself could have a typical monotonic dx^-y^ 
symmetry. 

The clue to understand the electron doped cuprate 
comes from two doping-dependent FS as revealed by 
ARPES [ll|, These are well explained in terms of 
the k-dependent band-folding effect associated with an 
AF order which splits the band into upper- and lower- 
branches [H, [III- In the SC state, the QPs could pair 
each other within the same band that leads naturally to 
a two-band/two-gap model. The two-gap model gives a 
unified explanation for the upward feature near Tc and 
the weak temperature dependence at low T in superfluid 



density ps [Ta|- It is also supported by Hall coefficient 
and magneto-resistance measurements [IgI [TtI [isj . 

Raman scattering has the potential to probe different 
regions of the FS. It has been shown by Lu and Wang 
that SC and AF orders in electron-doped cuprates are 
disentangled in Raman spectra. In our earlier calcula- 
tion on Raman spectra [l9| , we have also proved that the 
Raman shift for electron-doped cuprates is mainly de- 
termined by their pair breaking associated with different 
pieces of the FS. The AF order can cause a vertex correc- 
tion and enhance the spectral weight, but nevertheless, 
it does not change the Raman symmetry. In particu- 
lar, near the optimally-doped regime, the frequency of 
-^2o peak appears to be higher than that of Big peak 
It seems indicating that the SC gap deviates from 
the monotonic dx'^-y^-'waye. However, it has been shown 
that it is indeed two monotonic dr^2 _y2-w&Ye gaps, asso- 
ciated with a and /3-band FS respectively, which leads to 
a good description for it [19]. 

In this paper, spin dynamics is explored to further 
test the two-gap model for the electron-doped cuprates. 
Spin fluctuation is observable by inelastic neutron scat- 
tering (INS), and is confirmed to be intimately connected 
with the pairing mechanism in hole-doped cuprates. 
In single (CUO2) layer hole-doped cuprates such as 
La2_a;Sra;Cu04, the magnetic peak is always incommen- 
surate and their incommensurability is robust against the 
frequency change [20I W^. (Strong commensurate peak 
at momentum Q = (tt, tt) and some particular resonance 
frequency has been observed in multilayer YBa2Cu307 
and Bi2Sr2CaCu208 though [l^].) In current single- 
layer electron-doped cuprates such as Nd2-a;Ce2;Cu04 
(NCCO), in contrast, commensurate peak at Q is ob- 
served both in the SC and normal states [23l [2^. [25l. [26j . 
These commensurate peaks survive over a wide frequency 
range. It will be shown later that the commensurabil- 
ity of these magnetic peaks is a natural outcome of the 
band nesting, and their robustness is actually incorpo- 
rated into the existence of two separate bands. Spin 
dynamics has been theoretically examined in various as- 
pects for electron-doped cuprates lately [13, [H, [2^, [s^l • 

Based on a mechanism making use of the strong AF 
spin fluctuation, a pairing model will be proposed for the 
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electron-doped cuprates. Analogous to the hole-doped 
side, the SC gap function of the electron-doped cuprates 
is thus naturally to have the d^2_j^2 symmetry in the 
whole doping range. When the AF order is significant, 
it makes a big split between the two bands and conse- 
quently the nonmonotonic dx2_y2-wave like of the gap 
is satisfactorily explained. Of equal importantance, this 
model gives a unified picture for the Raman scattering, 
Ps{T), and INS in electron-doped cuprates. 



II. THE MODEL 

We start with a phenomenological superconducting 
Hamiltonian 

^ = E [^k/L/k,. + Ak (/,I,T/lka + /-k,i/k,T)" 

k,o- 

-2 Jm^'a(/^^,/k+Q,. + h.c.) - f, (1) 



originated from a t-t'-t"-J model. The slave-boson trans- 
formation and spin-density-wave mean-field approxima- 
tion are undertaken. Here /j[ ^ {fk,a) is the fermionic 
spinon creation (destruction) operator, m — (— l)*(S'f ) is 
the AF order, /z is the chemical potential, and 

£k ~ {2\t\d — J x) {cos + cos ky) — it' 6 cos kx cos ky 
-2t"S{cos 2kx + cos 2ky) (2) 

is the independent particle dispersion with S the doping 

concentration and x ~ ifjafj'^) ^^^'^ uniform bond order. 
The prime denotes that momentum summation is over 
the magnetic BZ only (— tt < k^ ± ky < tt) . The SC gap 
function is given self-consistently 



Ak = ^F(k,k')(/t,^/lk,,), 



(3) 



where V(k, k') is the pairing potential. Using the unitary 
transformation 



/k,(T 
fk+Q,a 



COS t^k cr sm f^k 
a sin 0k cos 



ak,<T 
/3k,. 



(4) 



with ^k being defined by tan 20k = 4Jm/[(T(£k+Q — Ek)], 
Hamiltonian ([1]) can be transformed as [3l|, [S^] 



H = 



E' [^ki^L^k^ + AkaJ^-r/^ki + '-kT^ki) -Mi, (5) 



ko-i 



"kT'-ki 

where I = a, P, fia + fJ-fj = /z, and 

Ek + Ek-I-Q 



ki 



^^/l£h±Q_£il)!+4J2^2 (6) 



corresponding to the QP (with only the AF order in it) 
dispersions of the two (a and /3) bands. 
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FIG. 1: (a) Separate a- and /3-band FSs of the electron-doped 
cuprates with an AF order, (b) Gap function Ak(0) with <j!> 
measured along the FS. Dash line: a nonmonotonic d^i_yi- 
wave Ak = Ao7k with 7k = sin 2(\> -\- a\ sin -I- a2 sin 6(/f> 
proposed in Ref. [33]. Here Aq = 33.5 cm~^, oi = 0.42, and 
a2 = 0.17. Solid lines: the piecewise monotonic dj.2_y2-wave 
Ak = Ai7k with 7k = sin 20 proposed by us. For < < 0/3, 
A; = A^ = 70 cm~\ while for (j>a < (p < 7r/4, A; = A^ = 25 
cm~^. The two gap amplitudes, Aq and Ap, are calculated 
using ([S} as a function of doping (c) and temperature (d). See 
Table I for parameters. 



Throughout this paper, |t| = 0.326 eV is taken as 
the energy unit together with t' = 0.3, t" = —0.2, and 
J = 0.3. Other doping-dependent parameters which 
agree well with the band FS are listed in Table HI 

Fig. [IK a) shows the FSs for a typical optimally-doped 
sample. The AF correlation splits the continuum FS 
into two pieces of sheet, in which a band is crossed by 
the Fermi level in the antinodal region, while f3 band is 
crossed by the Fermi level in the nodal region. Corre- 
spondingly small FS pockets appear around (±7r, 0) and 
(0, ±7r) for the a band, while separate FS pockets appear 
centered at (±7r/2, ±7r/2) for the (3 band. 



III. GAP SYMMETRY: MONOTONIC VS. 
NONMONOTONIC d^2„y2-WAVE 

As mentioned previously, ARPES ^] and Raman Q 
reveal the gap function Ak in electron-doped cuprate SC 
to have a non-monotonic (ia;2_j,2-wave like. For example 
in Ref. [ll], Ak = Ao7k with 7k = sin 20 -I- ai sin 40 -|- 
02 sin 60 (0 being the angle measured relative to the di- 
agonal on the FS) was used to simulate a non-monotonic 
dx2-y2-wave gap [see also Fig. [Ijb)]. In this kind of ap- 
proaches, parameters ai and 02 are doping dependent. 

As far as SC gap is concerned, it is physically more 
appealing that 7k remains the same for the entire doping 
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TABLE I; Parameters used to calculate and A/3. 
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m 


-X 




M/3 


0.15 


0.178 


0.15 


0.005 


0.078 


0.155 


0.169 


0.16 


0.007 


0.079 


0.16 


0.160 


0.17 


0.008 


0.079 


0.165 


0.150 


0.18 


0.009 


0.079 


0.17 


0.040 


0.20 


-0.040 


-0.036 


0.18 


0.010 


0.20 


-0.034 


-0.035 



range, so long as the pairing mechanism remains the same 
(no quantum criticahty occurs). Based on a mechanism 
induced by the strong AF spin fluctuation, we propose 
the foUowing piecewise model for the pairing potential 

{.9a 7k 7k' 7 for k, k' on a-band FS, 
5^7k7k' , for k, k' on /3-band FS, (7) 
0, otherwise, 

where 7k = sin2(/), and ga and gf3 are the two coupling 
constants. Inspired by the hole-doped side, it is promis- 
ingly to have 7k having the monotonic (1^2 _y2 symme- 
try. This is strongly supported by the INS experiment 
and will be elaborated later. Nevertheless, there are two 
dx2_y2-weiYe gaps, possibly with different amplitude, for 
the current electron-doped side. When AF order breaks 
down with increasing the doping, the two bands [Eq. ^] 
will eventually merge into a single one. In this regime, 
the behavior of the electron-doped cuprates is expected 
to be very similar to that of the hole-doped ones, with 
one single monotonic dx2_y2-wave gap Il9[. The latter is 
confirmed by the Raman experiment [6| . 

The piecewise feature in V^(k, k') comes naturally for 
the AF spin fluctuation induced mechanism (see Sec. IV). 
Only QPs within the same band favor the pairing asso- 
ciated with the Q wavevector . This is also supported 
by the superfluid data which unambiguously reveals that 
QPs in a (or /3) band FS pair each other to form SC 
QPs. But no SC QP forms in the region where FS is 
absent flE\. 

Substitution of ([7]) and ^ into ([3]), one obtains the 
self-consistent gap equation respectively for each band 

kei FS \ -D / 

where Ei^i = i/^^T+^ki" ^ki = Ai7k. In practice, 
the k sum in ^ can be effectively extended to the whole 
MBZ because contribution due to the k points distant 
from the corresponding Z-band FS is negligible. Fig. [TJb) 
shows an example of piecewise monotonic d^2_y2-wave 
Ak/, compared with a nomonotonic one. The region be- 
tween ((/)/3, (f)a) is where FS is absent and no SC gap as- 
sociated with. 

Shown in Fig. [2c) are the two gap amplitudes, 
and A^, calculated at T = with various doping levels. 
As the doping 5 decreases (and hence the AF order m 
increases), the ratio of A^/A^ increases along with the 



gapped ((/)/3, (/)„) region opens up. This manifests the non- 
monotonic d^2_j^2-wave like gap nicely. Both A^ and A^ 
decrease as doping increases. At over doping [8 > 0.17), 
m approaches zero and FSs join to one piece, Aq, and A^ 
match. These consistent results give strong support to 
the model pairing ([7]) . The parameters used are listed in 
Table m Coupling constants, ga and (7/3, are fixed at 0.34 
and 0.62 respectively, that give the best fit for optimal 
doping {5 = 0.15). Fig. [T]^d) displays the temperature 
dependence of A„ and A^j {5 = 0.15). The SC T^, de- 
termined by the higher of the onset temperatures that 
make Aq or A_a vanish, is found to be about 25 K. The 
two onset temperatures, differed by 4K or so, are in good 
agreement with the upward curvature observed in ps near 
Tc M- 

IV. DYNAMICAL SPIN SUSCEPTIBILITY 

The dynamical spin susceptibility, which comes from 
the particle-hole excitations, is given by 

X°(q,T) = l(T,5^(r)5^q(0))o, (9) 

where S^ = \ Y.^,^^ ^/k+q/2,fT/k-q/2,fT is the spin-density 
operator with cr the spin index. Using the transformation 
(jll) and Fourier transforming x°(qi ''') into the Matsubara 
frequency space, one obtains 

X°(q,iw„) = -^Xl ^^■M'xhi'i^^i^n). (10) 
k,ir 

Here i^k,(;'(q) = {1 + e/;/ cos[2(6lk - 6'k+q)]} with eu> = 
1 (-1) for I = l'{l ^ l') and 

xiwi^'i^n) ^ -■^^[tJ;(k,w„)C/;/(k-f q, zi/„ -hzw„) 

+ewTi{k, ivn)^l (k + q, w„ + iujn)] 

with Qi and J-i the single-particle normal and anomalous 
Green's function of band I. Considering the AF vertex 
correction under the random-phase approximation, one 
then has the renormalized spin susceptibility xi^-i i^n) — 
X°(q,«w„)/[l + J^J(q)x°(q,«^^«)], where J(q) = cos.{q^)^ 
cos{qy) and 1/ is the coupling strength. The INS intensity, 
/(q, cj), is proportional to Imx(q, it^^n ui + iO^). As a 
matter of fact, vertex correction leads to enhancement of 
the spectral intensity, but giving no qualitative change in 
the lineshape. 

Fig. [2ja) shows the calculation of /(q, lj) with differ- 
ent frequencies (w = 3, 5, 10 meV) at T = and optimal 
doping {6 — 0.15). The momentum is scanned along the 
direction of (0,0) (tt, tt). For easy comparison, the 
coupling strength 1/ is chosen to be 0.627, same as that 
used in Ref. [2g|. The most remarkable feature is that 
I{q,uj) is commensurate for a; < 5 meV, consistent with 
the INS measurements of Yamada et al. ^31 . In the case 
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FIG. 2: (a) Constant a;, q-dependent INS intensity, 7(q, oj), 
calculated for optimally-doped (& = 0.15, Tc = 18K) sample 
in the SC state (T — > 0). The smearings are taken to be 
= = 20 cm~^. (b) Energy contours of i5k = 2.5meV 
(5meV) shown in the second and fourth (first and third) quar- 
ters of the BZ. The double-arrow lines denote the correspond- 
ing nesting wavevectors. (c)&(d) Comparison between theo- 
retical calculations (solid lines) and experimental data (dots 
with error bars) of /(Q,t^) on SC NCCO (5 = 0.15). Tc = 18 
and 25 K for (c) & (d). 



of higher a; = 10 mcV, in contrast, /(q, w) becomes m- 
commensurate. The above theoretical results are in great 
contrast to those obtained based on a one-band model 
[28| . where spin response is found to be incommensurate 
at lower frequencies but shifted to be commensurate at 
higher frequencies. 

The switch from a commensurate to the incommensu- 
rate peaks upon frequency increase can be understood 
in terms of the band nesting effect. As illustrated in 
Fig. IHb), two sets of energy contours: = 2.5 meV 
and 5 meV are plotted respectively in the second and 
fourth and first and third quarters of the BZ. When en- 
ergy is low, only /3 band opens up a contour, and the flat 
(nesting) portion of the energy contours show a thin strip 
near the MBZ border. The corresponding wave vector 
q (double arrow) which brings the nesting portion into 
good alignment with its partner in the other quadrant 
equals to Q. Consequently, the nearly degenerate exci- 
tations give a commensurate peak. When energy is high, 
in contrast, both a and P bands open up a contour. In 
addition to Q nestings, the most contribution may come 
from the incommensurate Q ± 5 nestings [see Fig. (^fb)]. 

The exact or near (S is small) Q spin fluctuations may 



assist the QPs to form the d2.2_j,2-wave pairing. This 
is the scenario widely believed for the hole-doped side. 
Based on this pairing mechanism and taking possible Q 
connections into account [in view of Fig. 2(b)], the valid- 
ity of the model pairing potential ([7]) is justified. 

One can examine the intensity at q = Q more care- 
fully. At low w, /(Q,w) is weak, indicating that a 
spin gap opens up. /(Q,a;) will reach its maximum at 
Lu = 2A with a pairing-breaking gap estimated to be 



|A|- 



j^Tfii^Uo i-Kpi d(j) + J^^ [Akal When 
u! is higher, nesting portions move out of the MBZ 
boundary, and consequently I{Q,uj) starts to dimin- 
ish. A good agreement between the theoretical calcula- 
tion and experimental /(Q,w) is obtained and shown in 
Fig.[2i;c)&(d). In Fig.Etc) wi_th Tc = 18K, = 22cm-i 
and A/^ = 38cm^^, and 2A ~ 5.2 meV is obtained. 
While in Fig. ^d) with = 25K, A„ = 32cm- ^ and 
Ap = 55cm-^, and 2A ~ 7.6 meV is given. It is noted 
that the above A^ and A^ are taken exactly the same as 
those led to good fits for the Raman scattering [Toj . 

Recently, high-energy spin excitation of INS experi- 
ment is also reported [26] . The energy taken in those ex- 
periment is far above 2 A; 10 meV), beyond the scope 
of the present paper. In such case, the excitation leads to 
a spin- wave-like ring [2Q] and the effect of pairing break- 
ing is weak. A two-dimensional AF Heisenberg model 
including nearest (Ji), next-nearest (J2), and next-next- 
nearest (J3) couphngs should be used to interpret the 
experiments. 



V. SUMMARY 

In summary, the gap symmetry of electron-doped 
cuprate superconductors is studied based on a two-gap 
model. Considering a mechanism induced by the AF spin 
fluctuation, a piecewise pairing potential is proposed to 



account for the observed nonmonotonic d^. 



-wave fea- 



ture of the gap. Dynamical spin susceptibility is calcu- 
lated and shown to be in good agreement with the ex- 
periment. This gives a strong support to the proposed 
pairing model. 
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